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The qubit (or a system of two quantum dots) has become a standard paradigm for studying 
quantum information processes. Our focus is Decoherence due to interaction of the qubit with its 
environment, leading to noise. We consider quantum noise generated by a dissipative quantum bath. 
A detailed comparative study with the results for a classical noise source such as generated by a 
telegraph process, enables us to set limits on the applicability of this process vis a vis its quantum 
counterpart, as well as lend handle on the parameters that can be tuned for analyzing decoher- 
ence. Both Ohmic and non-Ohmic dissipations are treated and appropriate limits are analyzed for 
facilitating comparison with the telegraph process. 

PACS numbers: 03.65.Yz, 05.40.-a, 42.50.Lc 



I. INTRODUCTION 

Environment-induced decoherence is a vexed issue in 
confronting the challenges to quantum information pro- 
cesses. Simple model calculations thus assume signifi- 
cance, to grasp the underlying parameter-regimes that 
can be manipulated, in order to minimize the effect of the 
environment. In quantum optics and solid state physics, 
the device that has gained popularity in recent times is 
a qubit (which can be approximately represented by a 
pair of two quantum dots)^''^. This system is equivalent 
to a superconducting Josephson junction^ or a spin-1/2 
NMR nucleus^. A qubit can be described by a two-level 
Hamiltonian, the quantum mechanics of which is rather 
straightforward. Environmental influences can also be 
easily incorporated in this Hamiltonian. There have been 
two distinctive attempts in modeling the environment, ei- 
ther in terms of (i) a classical stochastic process, such as 
a Gaussian or a telegraph one, or (ii) an explicit quan- 
tum collection of bosonic oscillators. One of the main 
objectives in this paper is to seek a unification of these 
two apparently disparate approaches. 

We consider an electron hopping between two single 
level quantum dots, coupled via a tunneling coefficient A 
and an asymmetric bias e. Denoting by \L) and \R) the 
left and right dot states, the qubit Hamiltonian, indicated 
by the subscript q, can be written as, 

n, - A(|L)(i?| + \R){L\) + ei\L)(L\ - \R){R\) (1) 

If we introduce the so-called bonding and antibonding 
states as the spin 'up' and spin 'down' states of the Pauli 
matrix cr^, denoted by |-|-) and |— ) respectively, as 

\+) = ^i\L) + \R)), 

h)^^i\L)-\R)). (2) 
Hq can be cast into more familiar notation: 



'Hq = Aaz+€ax, (3) 

where is the other Pauli matrix that is off-diagonal in 
the representation in which az is diagonal. 

In order to incorporate the effect of the environment 
on the qubit our strategy is to expand the Hilbert space 
to re-express the system Hamiltonian as 

Hs=Hq+TziCA^z+Ce<Jx), (4) 

where Tz is yet another Pauli matrix, which could repre- 
sent a spin-1/2 impurity, for instance. The idea is to cou- 
ple the T^-system to a bath of bosonic oscillators, thereby 
causing fluctuations in . The resulting full Hamiltonian 
can be written as 

H^H^+T^Y. 9k{bk + bl) + Y^ iOkblbk. (5) 

k k 

Here 6[.(6fe) are boson creation (annihilation) operators, 
g'l^ are coupling constants and ujk is the harmonic oscil- 
lator frequency of the /c*'' mode. 

It is interesting to physically assess the effect of the 
coupling term (proportional to gk) on the system Hamil- 
tonian Us- Because is purely off-diagonal in the 
Tz- representation, it would cause 'spin-flips' in be- 
tween two allowed values -1-1 and -1, as in the celebrated 
Glauber model of Ising kinetics^. These flips would oc- 
cur with an amplitude fleld that would be proportional 
to gk and time- varying bosonic operators bl,{t){bk{t)), in 
the interaction representation of the last term in Eq. (5). 
The net effect on Jis is a 'quantum noise', encapsulated 
by Tz{t). In a suitable limit, when the latter, i.e. Tz{t), 
could be replaced by a classical noise r]{t) that jumps at 
random between ±1, the system Hamiltonian would be 
stochastic, given by, 

Hs{t)=Hq+v{t){CA<yz+Ce(Tx), (6) 
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where r]{t) is a two-state jump process or a Telegraph 
Process. While we will present below a fully quantum 
mechanical treatment of Eq. (5), our aim will also be 
to set conditions under which a classical stochastic de- 
scription via a telegraph process of ri{t), as expounded 
in detail by Tokura and Itakura^, would ensue from the 
quantum formulation. 

The form of the Hamiltonian in Eq. (5) falls under 
the general scheme of a system-plus-bath approach to 
nonequilibrium statistical mechanics, in which the full 
Hamiltonian is written as, 



n = ns + ni + 'HB, 



(7) 



where Hi is the interaction Hamiltonian and 'Hb is the 
bath Hamiltonian. In the present instance, 



Hi = T^'^gkih + bl), 

k 

Hb = ^Wfc&i^fc- 



(8) 



the much-studied Gaussian stochastic processes. Simi- 
lar spin-boson Hamiltonians, akin to Eq. (5), have been 
considered in the past^^, though the comparison has not 
been dealt with, to the best of our knowledge. 

With these preliminaries the paper is section wise or- 
ganised as follows. In Sec. II we outline the mathemati- 
cal steps for a resolvent expansion of the averaged time- 
development operator in the Laplace transform space, 
the average being carried out over the Hilbert spaces of 
and "Hfi. The latter is expressed in terms of a 'self-energy' 
whose form is provided. In Sec. Ill explicit results 
are presented for the much explored Ohmic-dissipation 
model and comparisons are drawn with the classical tele- 
graph process, results of which are given in the Appendix. 
For facilitating comparison we discuss fluctuation in bias 
and hopping separately. In a subsection of Sec. HI, we 
indicate results for non-Ohmic situations as well. Finally 
in Sec. IV, we discuss the issue of partial decoherence, 
recently studied by us in the context of the telegraph 
process"'^^, to put it under the perspective of quantum 
noise. The section V contains a few concluding remarks. 



As wc will be employing a Liouvillean approach to the 
dynamics governed by "H, it is useful to explain the no- 
tation. The Liouville operator £, associated with %, is 
defined by 



CA=-[n,A], 



(9) 



where A is an arbitrary but ordinary operator. Thus a 
Liouville operator yields an ordinary operator when it 
operates on an operator (as in the right hand side of 
Eq. (9)), just as an operator, operating on a wavefunc- 
tion, yields another wavefunction. We will designate the 
'states' of a Liouville operator by roimd brackets: |), in 
analogy with the Dirac ket vectors: |) for the wavefunc- 
tions. The decomposition in Eq. (7) further allows £ to 
be split as 



C = Cs+Ci+Ci 



(10) 



The Liouville dynamics effected by C on the density 
operator p, leading to a non-Markovian master equa- 
tion, in the context of the decoherence of a qubit (as 
well as a collection of qubits) have been extensively stud- 
ied by Kurizki and collaborators^. They have also con- 
sidered the additional influence of external time-varying 
fields for controlling decoherence, dynamically. While 
we will be employing a very similar approach, our treat- 
ment will differ from Gordon etal®. The latter utilize the 
Born approximation, valid for weak coupling constants 
Qk, whereas we will consider strong coupling in which gk 
will be treated to all orders. 

Such strong coupling considerations are the hallmark 
of functional integral treatments of quantum dissipative 
systems eg. that of a spin-boson Hamiltonian^". Fur- 
thermore our focus in the classical limit of the quan- 
tum noise will be a telegraph process as opposed to 



II. THEORETICAL FORMULATION 

In order to incorporate strong-coupling effects it is con- 
venient to transform the Hamiltonian T-L in Eq. (5) with 
the aid of a unitary transformation defined by the oper- 
ator: 



exp 



exp 



TT 



The transformed Hamiltonian is given by 

n = sns-\ 



(11) 



(12) 



What the second term in S docs is to cause a rotation in 
the T-space by an angle of 7r/2 about the y-axis, in the 
anti-clockwise direction such that Tz and — r^. 
The first term in S then eliminates the interaction term 
but puts the onus of coupling on l-Ls itself (besides 
generating an innocuous 'counter term' that is constant, 
and can be dropped). The net result is 



where 



B± = exp 



(13) 



(14) 



The transformed Hamiltonian 71 is now endowed with a 
new interaction term that may be written as 

1 



■Hi 



(t-B+ + T+B_) (CaCT. + Ce<7x) . 



(15) 



3 



The important point however is that any perturbation 
treatment of Jij is tantamount to treating the coupHng 
to all orders as Qq now occurs in the exponent as is evident 
from the structure of B± (cf. Eq. (14)). The equation 
(13) is again of the generic form of Eq. (7) except that 
the 'system' is now the qubit itself as Hs is replaced by 

With reference to the Liovillean in Eq. (10) the equa- 
tion of motion for the density operator p reads 



-iCp{t). 



(16) 



Unlike other approaches^^^ we find it convenient to work 
with the Laplace transforms, defined by 



From Eq. (16) then 



Jo 



dte-''^p{t). 



(17) 



piz) = {U{z))p{t = 0), 
U{z) = {z + i£)~\ 



(18) 



We employ the usual factorization approximation that at 
t — the qubit and the bath are decoupled^, so that 



p{t = 0) = Pg{0)®pB, 



(19) 



where Pq{0) is the qubit density operator at time i = 0, 
Pb = gxp(—I3'Hb)/Zb, is the canonical density opera- 
tor for the bath and Zb is the corresponding partition 
function. Underlying the prescription in Eq. (19) is the 
assumption that the bath always remains in quilibrium 
at a fixed temperature T while the qubit evolves from an 
arbitrary state. For most of our derived expressions we 
shall assume that the electron is on the left dot at f = 



I.e. 



p,(0) = \L){L\. 



(20) 



Translated to the 'bonding' and 'anti-bonding' basis, this 
implies 



P,(0) = ^(l + a,). 



(21) 



(We shall return in Sec. IV to a more general initial con- 
dition.) The Laplace transform of the so-called reduced 
density operator can be written as 

pn{z) = Trr Tr^ [{U{z)) PB]Pg(0), (22) 



where Trr{.-.) denotes trace over the eigenstates of Tz 
whereas TrB{---) specifies trace over the bath (i.e. eigen- 
states of Hb)- Because the trace is invariant under the 
unitary transformations (eg. in Eq. (11)), Eq. (22) can 
be equivalently expressed as 

pniz) = Trr®TrB [{u{z)) Pi3]p,(0). (23) 
where the tilde on U{z) implies that the corresponding 
time-evolution is now governed by U of Eq. (13). 

Denoting the Liouville operator associated with i-Lj as 
£/, we have upto second order in 



Pr{z) 



z + i. 



(24) 



jCiPb 



where the self-energy '^(z) is 

J2{z) = Trr Ttb [Ci {z + i{Cg + Cb)) 

]25) 

The z-dependence (or frequency dependence) of X](z) 
implies that the underlying dynamics is non-Markovian, 
in general. We emphasize once more that though the 
self-energy is computed to second order in £/, the the- 
ory is valid to arbitrary orders in the original coupling 
constants gk- Hence, the resultant master equation in 
the time domain that emanates from Eq. (24) is of 
more general validity than the one obtained in the Born- 
approximation^ . 

Our next step is to write the matrix elements of the 
self-energy '^{z). It is clear that ^^^(z) is a Liouvillean 
in the cr-subspacc alone because the t and the bath vari- 
ables are averaged over, in Eq. (25). Hence, '^{z) has a 
4x4 matrix representation in the 'eigen basis' \pv) (note 
the round brackets used for 'states' of the Liouvillean), 
wherein |^), \v),... are the eigenkets of cr^. Clearly then, 
the rows and columns of the 4x4 matrix would be la- 
belled by I -|- -h), I ), I H — ) and | — h), respectively. 

Additionally, while performing the traces in Eq. (25), we 
would need to sum over the eigenstates of and 'Hb, s& 
in the following: 



Tz\Sl) = Si|si), 

HbW) = En^\nb), 



(26) 



etc, With these explanatory remarks about the notation, 
the matrix elements of J2{z) turn out to be generalized 
forms of those given in Ref. [13]: 
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^ {v'nhSi\Hi \r1nls2) {'nnls2\ Hi \vnbs{) 

r, Z + i{E^-Er,)+i[Er,,-E„ij 

{lJ,nbSi\Hi |/u'n^S2) (z^'n^S2| -ff/ Wmsi) 

z + i {E^> -E^) + i (i;„i - 
{linbSx\Hi \^i'nls2) {v'nls2\Hi \vni,si) 

z + i{E^-E,,) + i(E^^-E^.^ 



]{nb\pB\nb). 



(27) 



From the structure of T-Li in Eq. (15) it is evident that the bath variables enter only through the operators B±. 
Therefore, when wc sum over the bath indices Hb, n^, etc., with the Boltzmann weight (n(,|ps|n(,), we end up with 
correlation fucntions of B±, as in the following: 



^ij{t) = Ttb {pBBdO)Bj{t)) = {B,{0)Bj{t)) = +or-). 

As it happens, the only non-zero components of the correlation function are 



(28) 



$+_(*) = $_+(t) = (i3± (0)5^(0) 



exp 




coth 



( 1 — cos ojgt) +i sin ojqt 



(29) 



An essential attribute of quantum dissipative system is that the system He is taken to the limit of an infinitely 
large number of bosonic modes. This implies that the discrete sum over q has to be replaced by an integral over 
continously varying frequencies oj with an appropriate weight called the 'spectral density' J(w), thus 



</)(r)=exp^-2 £ d^"^ 



where 



coth ( ^ ) (1 — coswr) + isinwT 



(30) 



(31) 



r 



III. EXPLICIT RESULTS AND COMPARISON 
WITH TELEGRAPH NOISE 

A. Ohmic Dissipation 

In the literature on dissipative quantum systems one 
model for the spectral density J{uj) that has received 
the maximum amount of attention is what gives rise to 
Ohmic dissipation^*. In this, J(w) is assumed to have a 
linear frequency dependence with an exponential cut-off 



J(w) = Kuj exp ( — 



(32) 



where K \s& phenomenological damping parameter that 
measures the strength of the coupling with the heat 
bath. This form allows for an analytic expression for the 
Laplace transform as defined in Eq. (17), for j3wc >> 1 
as 



27r 



(2K-1) 



T{1-2K)T{K + zI5/2t:) 
Til-K + zp/iir) 



4>{z) 

UJc \ l->^c J L {1 — JTL -f ^IJ/ ■i'/i j 

(33) 

where r(..) denotes gamma functions. 

The significance of the condition /3ojc >> 1 can be as- 
certained from the fact that the quantity Htq = /3 (Note: 
h has been set to unity, hitherto) sets the time scale over 
which quantum coherence is maintained. Thus, for all 
frequencies greater than the cut-off uic, quantum fiuctu- 
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ations remain strictly coherent. Before we take up the 
comparison with the telegraph noise it is pertinent to 
point out that it is meaningful to think of a quantum 
bath as classical only in the incoherent regime i.e. over 
time scales larger than Tq. In other words, the classical 
limit of the bath obtains for frequencies uj « l/r^, i.e. 
the temperature is significantly large, but not so large as 
to violate the condition ujcTq » 1, if we want to restrict 
our discussion within the domain of validity of the ana- 
lytical result of Eq. (33). If the restriction wTq << 1 does 
apply, the frequency or z-dependence in the arguments of 
the Gamma functions in Eq. (33) can be dropped yield- 
ing 



A = 0(z = 0) = 



27r 



{2K-1) 



r(l - 2K)T{K) 



(34) 

The parameter A will turn out later to be related to the 
jump rate of the telegraph process. 

With this background we are ready to provide explicit 
results for the different elements of the density matrix in 
order to examine decoherence as well as to draw compar- 
ison with the results for the telegraph noise. Because the 
latter have been extensively studied recently by Itakura 
and Tokura'', we will follow their lead in ignoring the 
bias term in the original qubit Hamiltonian Hq, i.e. set 
e = though it should be evidently clear that the formal- 
ism can easily embrace asymmetric cases (e 7^ 0) as well. 
With this, the relevant Hamiltonian for further consider- 
ations can then be written as (cf., Eqs. (3) and (13)) 



where 

ai(z) = 2C[4>{z+) + 4>'{z-)],a2{z)^2CnH^-) + i'{z+)] 
a^iz) = 2(:^[^{z) + 4>'{z)], z±=z±2zA. (38) 

Here (j){z) and (p'iz) are the Laplace transforms of the 
correlation functions (f>{t) and 4>{—t), respectively. 

The matrix element of (^z + iCq + X](z)^ , required in 

Eq. (24), is obtained from a combination of Eq. (27) 
and the definition of the matrix elements of the Liouville 
operator £q a la Eq. (A. 7). The resultant matrix is of a 
block-diagonal form and can be easily inverted to yield 

{PTz{z)) + + = {pTz{z))__ 



2z' 

z + 2iA + 442 



[z) + 0'(z) 



z z-f 4(2 U(z) +0'(z) +4A2 




FIG. 1: Im[{pn)LR (i)] is plotted for ^ = 100, k ^ 0.45 
and A = 2 



n = Acr,-- 



k 



k"k 



(35) 



For ease of discussion and for remaining close to the 
treatment of Itakura and Tokura we take up the cases 
of fiuctuation in bias (Ca = 0) a-nd fluctuation in hop- 
ping (^e — 0) separately below. 



B. Fluctuation in Bias ((Ja = 0) 



Villi 
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The appropriate interaction Hamiltonian that has to 
be substituted in th expression for the self energy '^{z) 
in Eq. (25) is now given by 



Hi ^--C,a^{T-B+ + T+B. 



(36) 



Carrying out the summations implied in Eq. (25) we 
arrive at 



ai{z) — ai(z) 

-02(2) 02(2:) 

03(2) -03(2) 

-03(2) 03(2) 



, (37) 



FIG. 2: {pn)LL (0 is plotted for ^ ^ 100, k = 0.45 and 

A = 2 

It is of course possible to go back to the dot-basis with 
the aid of Eq. (2). Results are shown in Figs. 1, 2 and 
3 in which Im[pLii{t)], pLL{i) and pRR(t) are plotted 
versus time t for a certain choice of parameters. While 
ImlPLRit)] oscillates and decays to zero as i — > 00, the 
diagonal elements, proportional to the populations of the 
two dot states, oscillate and settle to the common value 
of 0.5, as we expect for full decoherence. If we compare 
Figs 1-3 with Fig A.l for the telegraph process below, it 
is evident that the quantum case exhibits persistence of 
oscillations over longer time scales. 
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FIG. 3: {pn)RR (0 is plotted for P = 100, k = 0.45 and 

A = 2 

C. Fluctuation in Hopping (i^e — 0) 




-0.5 



FIG. 4: {ptz)ll (*) (solid), {ptz)rr (0 (dashed) and 
Im {ptz)rr (t) (dotted) are plotted for P = 100, 
k = 0.45 and A = 2 



The appropriate interaction Hamiltonian that has to 
be substituted in the expression for the self energy 
in Eq. (25) is now given by 



Hi = {r'B+ + t+BJ) (40) 

Following the formalism developed in Sec. 11 and using 
Eq. (25), the self energy matrix '^{z) is given by 



E(-)- 







63(2) 

hi{z) 



(41) 



where 



h{z) 



4Ci[0'(z+) + (/-(z+)], 



b^{z) = 4Ci[0'(2-) + 0(z_)], ^± =z±2iA, (42) 

where (j){z) and 4>'{z) have been defined in the earlier 
section. ^ 

Using the form U{z) — (^z + iCs + ■• i* is evi- 

dent from the diagonal form of the matrix that the pop- 
ulation of the bonding and anti bonding state remains 
constant for all t > Q. Using U{z) the matrix elements 
are given by, 

{p-r)++ {z) 

(PTZ), (z) 



(pr)- (2) 

1 r 1 

2 



1 

2^' 



2iA + 6.S 



(43) 



However in the dot basis we expect the usual saturation 
to the common value of 0.5 for both left and right dot 
populations. Fig. (4) shows the variation of left and right 
dot populations with time and also the loss of coherence 
over time from the decay of Im{pLR{t)). 

We are now in position to make a connection between 
the results of quantum and classical regimes considered 
separately in Sec. Ill and Appendix A. There are two crit- 
ical distinctions between the two regimes. At the outset. 



the quantum behavior is non-Markovian, reflected in the 
z (or frequcncy)-dependence of the self-energy '^{z) in 
Eq. (25). On the other hand, in the classical stochas- 
tic case, the Markovian assumption is invoked, to be- 
gin with, as seen in Eq. (A. 2) below. The second 
crucial ingredient in the quantum situation, even after 
the Markovian approximation is incorporated, is in the 
temperature-dependence of the relaxation rate A. While 
the latter, in the classical domain, is usually endowed 
with an exponential dependence on temperature, with a 
barrier-activation in mind, it has a richer structure in the 
quantum regime, as discussed below. 

We mention below in Appendix. A. 2 that in the long 
time region the off-diagonal element of the system density 
operator follows an exponential decay with time constant 
T2 = I ^ From Eq. (43), as it can be clearly 

seen from the Laplace transform of p^ , the decay is 

exponential only if we assume 63(2) to be independent of 
z, wherein the time constant is just the inverse of i?e[63]. 
Thus the comparison with the classical telegraph process 
is meaningful only if one takes the Markovian limit of 
heat bath induced relaxation, at the outset. 

At sufficiently high temperatures ^iven by ksT » 
(a/K we can see from Eq. (33) that (f){z) can be approx- 



imated as 



2K-1 



r{i-2K)r{K) 
r(i--ff) 



which 



leads to a rather simple expression for 63 (for z = 0) 

in Eq. (42) as 63 = ^ cos (.X)(|^)^^-^ ^^Va-^T^ 
that is real. Hence A, which appears as the jump rate 
in the case of telegraph noise, discussed in the Appendix 
in the sequel, can be compared with the relaxation rate 
_Re[63] in the case of quantum noise, is now a function of 
K and /3. Also, as it is obvious from the expression of the 
relaxation rate that the dependence on temperature is 
not an exponential, rather it shows a power law behavior. 
However, it is pertinent to emphasize that these conclu- 
sions hinge on our assumption that quantum dissipation 
is Ohmic. In the more general case, discussed below in 
Sec. HI D, we will see that the temperature-dependence 
of the relaxation rate can in fact be exponential, under 
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certain special situations. 

At low temperatures and weak damping i.e. 'small' 



K, 4>{z) is given by (f){z) = 



z+27r//3 



where /z is a renor- 



malized form of tunneling frequency given by /i = 
4Ci(|^)2-^-ir(l - 2i^). Using this form of ${z), and 
Eq. (43) we can deduce, 



{pn)+- {z) 



r 4- 9i A _ 4iAcos(7rif)^^ 

^ -r ZtiA CO^/«12_i_/f A2 T 



4:ti^(TT /fi) COS (tvK) 



(27r//3)2+4A2 ^ (27r/;g)2+4A2 

(44) 

Prom Eq. (44), it is evident that the relaxation rate 

becomes ^^(2^/^)2^4/1^^^ ' which is far from an exponen- 
tial in temperature. Thus it may be stressed that even in 
the Markovian limit the temperature dependence of both 
the relaxation rate as well as the Rabi frequency is much 
more complex than an exponential, as is usually assumed 
for classical activated processes. 

These conclusions are similar to those obtained earlier 
in the context of spin relaxation of a muon, tunneling in 



a double well^^. 



D. Non-Ohmic case 

We had in the text presented results based on the 
-issumption of Ohmic dissipation i.e. choosing J{ix}) = 
^wexp (— ^)- The Ohmic model is a rather limited one, 
: lot applicable to phonons but is relevant to a bath con- 
^ listing of electron- hole excitations near the Fermi surface, 
as in metals. -'^■^'^^ In this subsection we will allow for a 
general spectral density function, and will be analysing 
the relaxation rate in the case of quantum noise. We 
will try to make a connection with the classical Tele- 
graph process and show that in the high-temperature 
limit, the relaxation rate indeed follows an exponential 
dependence on temperature with however quantum cor- 
rections, which was not the case in the Ohmic limit. Our 
starting point is Eq.(30) which upon clubbing the oscil- 
lating terms together, can be rewritten as, 



0(r) = exp 



dui 



J{co) 



coth 



l3u}\ cos(w(T-i/3/2)) 



sinh(if) 

By making a change of variable from r — i/3/2 to t in Eq. (45), ^{z) can be written as. 



(45) 



0(2) = exp 



die ^*exp 



-2 / duj 



J{u:) 



coth 



(3uj \ cos {ujt) 



2 J sinh(^) 



(46) 



In the high temperature regime^^, the integrand in the second term in Eq. (46), can be replaced by its short time 
limit as, 



^(.) = exp(-^)/_^^die-expL2^ 



Jiu;) 



tanh 



/3w 



+ 



2+2 



2sinh(^) 



(47) 



Also in the high temperature regime i.e. low P, tanh(a;) and sinh(a;) can be replaced by the argument i.e. x, and 
hence the integrations in Eq. (47) can be carried out easily to yield. 



(j){z) = exp (-/3(a + iz/2)) 



/OO 



e ^ , 



(48) 



where a = ^ duj^^^. Also (t>'{z) is obtained by replacing i by —i in Eq. (48). Thus the relaxation rate which is 
given by Re(53), for the quantum case in Sec. III.C, is obtained by imposing the Markovian limit i.e. putting z = Q, 
leading to 



63 = 4a(exp[-/3(a-A-^) 



/ rfwexp (-Aav^ip) + exp f-/3(a A - —) \ [ duexp (-Aau^/P) | . 

. m 



Finally the lower limit in Eq. (49) can be set to zero and follows a power law. Thus the dominant behavior of the 
thus the integral is just a function of temperature which relaxation rate is an exponential dependence on tempera- 



8 



ture, which is what we expect at high temperatures. This 
analysis therefore provides a regime where we can more 
meaningfully compare the relaxation rate A for the case of 
classical telegraph process with the Laplace transform of 
the bath correlation function in the quantum case. Note 
that the present treatment does not depend on any spe- 
cific assumption for the spectral density J{uj) but is more 
generally couched, irrespective of whether the phonons 
are acoustic or optic. 




IV. THE ISSUE OF PARTIAL DECOHERENCE 



FIG. 5: Re[{pTi)_^_ {t)] for quantum (solid) and 
classical (Dashed) cases are plotted versus time t. 



In the context of fluctuation in hopping, it is evident 
that the system Hamiltonian commutes with the inter- 
action Hamiltonian, which means there is no energy ex- 
change between the system and the bath. Such cases are 
important for studying partial decoherence* . Thus we 
generalize our discussion to a more general initial state 
of a qubit (than of Eq. (20)) given by, 



I*) 



e^^ sin I - ]\R). 



(50) 



In such a scenario it can be established that the den- 
sity matrix does not evolve to a completely mixed state, 
rather it approches a limiting value, which does contain 
the off-diagonal components. In this situtation coherence 
is not completely lost and the off-diagonal terms give in- 
formation about the initial state of the system. We note 
that the initial qubit information can be retrieved, ei- 
ther from experiments related to persistent current or by 
measuring the population of states. The former path is 
not applicable in our case as we have not allowed for the 
existence of an Aharonov-Bohm flux* . 

The density matrix (in \L) and \R) basis) at time t — 
is given from Eq. (50) by. 



Pk(0) 



cos2(6'/2) 
cos(6i/2)sin(6l/2)e*'^ 



cos(6'/2) sin(6'/2)e-^'^ 



sm 



'(0/2) 



(51) 

which can be translated to the bonding-antibonding ba- 
sis, via the transformation given by Eq. (2). The off- 
diagonal term {pti)lr{^ ~^ which measures coher- 
ence, does not go to zero unlike in the cases discussed 
in Sec. (III. C and A. 2). These particular symmetric cases 
of e = and Ce = show partial coherence which is im- 
portant for decoherence-free quantum computation pro- 
tocols. 

The above formalism is general and is valid for any kind 
of environment, as long as the commutator of the system 
Hamiltonian and the interaction Hamiltonian is zero. As 
shown in the cases of classical and quantum noises (Figs. 

5 and 6), the imaginary and real parts of [pti)^ {t) — > 

as t — >■ (X), whereas the populations of bonding and 
anti-bonding states stay constant, as is evident from Eq. 
(43). Thus, in general, the off-diagonal elements of the 
density matrix (in \L) and \R) basis) approach limiting 



k4\ A A A [\ f\ f. 



FIG. 6: Im[{pTz)^ (t)] for quantum (solid) and 

classical (Dashed) cases are plotted versus time t. 



values, emphasizing that the asymptotic state is not a 
fully mixed state: 



oo 



\ Re{{pn)LR (0)) 



(52) 

At this stage it is also worthwhile to underscore the 
importance of time scales involved in attaining partial 
decoherence. For this we define a function C(t) as: 



IM(.Pk)lr (oo)) - I'm{{pn)LR (0)) ' 



(53) 



This function is 1 at time t = Q and at time t = oo. 
When integrated over all times, we get a 'relaxation' time 
for substenance of decoherence as: 



f 

Jo 



dtC(t). 



(54) 



We carry out an explicit calculation for r in both classical 
and quantum regimes. Using Eq. (A. 14) and (43) and the 
initial density matrix given by Eq. (45), we can calculate 
iPT^)LR (*) to obtain: 



^h+l2lmiipn)LRm 
Im{{pn)LRm 



(55) 
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where, for the quantum case, 
h = 



(2A)2+&r 

(2A)2 + 6i' 



(56) 



On the other hand, for the classical telegraph process, 
(2A) (-A + k) (2A) (A + k) 



h 
h 



(2A)2 + (A±fe)2 (2A)2 + (A^)" 
(A±^)(_A + fc) (A_fc)(_A + fc) 



+ 



(2A)2 + (A±fe)^ (2A)2 + 
k = ^(A2-16Ci). 



A-fc 
2 



. 2 ' 



(57) 



where 63 is given by the z = Q limit of Eq. (42). It 
is evident from Fig. (6) that the partial decoherence is 
attained more rapidly in the case of classical noise. 



V. CONCLUDING REMARKS 

Our main emphasis in this paper has been a rela- 
tive assessment of quantum and classical noise sources 
attached to a pair of quantum dots or a qubit. The 
analysis has been made with the aid of a unified for- 
malism, comprising a rcsolvant expansion of an averaged 
time-development operator. The quantum formalism en- 
abled us to provide a microscopic meaning and detailed 
temperature-dependence to the phenomenologically in- 
troduced parameter of the relaxation rate, that appears 
in the classical case of a telegraphic noise. Our treat- 
ment of the quantum noise has included the much studied 
ohmic model of dissipation that characterizes electron- 
hole excitations off the Fermi surface, in a metal, as well 
as non-ohmic dissipation which covers both acoustic and 
optic phonons^^. It is eventually in the phonon model of 



dissipation that the usually assumed exponential temper- 
ature dependence of the relaxation rate (of a Telegraph 
process) is realized, making the comparison between the 
classical and quantum cases more direct. 

In the last subsection of the paper (see IV) we focussed 
on the important issue of partial decoherence that can be 
utilized for quantum computation, which has received re- 
cent attention^. This case was analyzed when fluctuation 
is ascribed only to the hopping term, yielding a situation 
in which the system Hamiltonian commutes with the cou- 
pling to the heat bath. Here, the comparison between 
classical and quantum noises is quite striking coherence 
persists over longer time scales for the quantum case. 
This attribute can be effectively exploited in the context 
of quantum computation, in which it is essential to be 
able to retrieve information about the initial quantum 
state, notwithstanding heat bath-induced effects. 
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Appendix A: Results for Telegraph Process 

1. Fluctutaion in Bias (Ca = 0) under Telegraph 
Process 

The stochastic Hamiltonian appropriate to Eqs. (35) 
and (36) is given by (cf., Eq.(6)) 



(Al) 



where ri{t) is a telegraph process. Concomitantly Eq. 
(24) is replaced by^^ 



Pn{z) = 



a,b \ j=± 



P«(0), 



(A2) 



where the 'stochastic states' |a), |5)... are associated with 
the two possible values ±1 of ri{t) (over which the sum- 
mations in Eq. (A. 2) are performed), pb is the a- priori 
probability of the occurence of the state |6), Cq is the Li- 
ouvillc operator accompanying Hq = Acr^, and Cj is the 
Liouville operator associated with an ordinary operator 
Vj (j=-|-l or -1) such that 



V± = ±QcT,. (A3) 

It remains then to define the 'stochastic matrices' Fj and 
W. The matrix Fj is a projection governed by 

{a\Fj\b) = SabSaj, (A4) 

whereas W is a jump matrix whose element (a|iy|&), for 
instance, is the rate of jump of the process r]{t) from 
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the state \b) to the state \a). The underlying Markovian 
assumption is manifest in the frequency independence of 
W. 

The properties of the telegraph process allow us to cal- 
culate the average over the stochastic indices a,b,... in 
Eq. (A. 2) in closed form, leading to 



Pn{z) = 



Uo{z + X) 



P,(0), 



1 - A[/o(z + A)_ 
where 

(Uoiz + A)) = ^pj- (z + A + tCg + iLjy^ 



(A5) 



(A6) 



A being the jump rate. 

In calculating the matrix elements of Uq{z + A) we ned 
to use the properties of Liouville operators given by^^ 

{^Jiv\^^^^i'v') = {^l\v^^i')5,,. ~ {v'\v,\v)6^^,. (A7) 

Using (A. 7), it is straightforward to obtain 



[z + iCq + iC±) 



where 



z 




±<. 


■ 




z + 2iA 





±iQ 







z - 2iA 











z 



z — z + X 



(A8) 
(A9) 



The matrix of Uo{z) is therefore of a simplified block- 
diagonal form from which matrix under the square brack- 
ets in Eq. (A. 5) can be evaluated easily. Combined 
with the initial condition in Eq. (21) we find, after some 
straight forward algebra, 



{pn)++ (z) 
(P7^)__ (z) 
and 

{Pn)+- (z) 
where 



2(q + b)- Xja^ - b^) 
2{a + b) - X{a^ - b^) 



2{c + b) ~ X(dd* ~ b"^) 



2(z2 + 2(2 + 4A2 



(.PTz)++ (0) 



iPn)+. (0) 

(All) 



D2 



z[z^ 




+ A2)] 








z[z^ 


+ 4(C2 


+ A2)] 




+ z{z- 


- 2iA)] 


z[z^ 


+ 4(C2 


+ A2)] 


(2- 


Aa)2- 


(A6)2, 


(2- 


Ac)(2- 


-Ac*)- 



FIG. 1: {ptz)ll (*) i^^M), {pn)RR {t) (Dashed), and 
Im[{pTi) {t}] (Dotted) are plotted versus time t for 
A = 1, Ce = 1-5 and A = 0.5 



With the chosen initial condition in Eq. (21), Eq. (A. 10) 

simplifies to {pb.{z))^^ = {pr{z)) = ^. It is of course 

a matter of ease to go from Eqs. (A. 10) and (A. 11) to the 
dot basis, with the aid of the transformation in Eq. (2) 
to rederive the elements of the density operator, results 
of which are presented in Fig. (A.l). The latter are in 
complete conformity with those of Ref. [7] and display 
full decoherence at infinite times i.e. {ptz)ll ^ iPTi-) 



and {pn)LR = as expected. 



Fluctuation in Hopping ((^j 
Process 



IRR 



0) under Telegraph 



The relevant stochastic Hamiltonian is now given by 
(cf., Eq. (6)) 



■Hs{t) 



(A13) 



In this case the fluctuating part of the Hamiltonian com- 
mutes with the system part of the Hamiltonian, a case 
which has received attention recently^2_ Jn this case now 
there is no energy transfer between the system and the 
bath, allowing for partial decoherence (see Sec. IV). 
In this case the calculations are simpler than in A.I as 
the matrix of Uo{z) is diagonal in the representation. 
Therefore, the diagonal elememts of pr{z) do not evolve 
at all from their initial values of one-half whereas the 
off-diagonal element becomes 



{pn{z))+_ 



2iA 



4Ci 



z-t- 2iA 



(A14) 



(A12) 



The time dependence of the off-diagonal term 

(p7j(z))_i (t) is obtained from the inverse Laplace 

transform of {pti{z))^ {z) to yield, for A > 4(a, 
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4yp 



i6ci 



A2 - 16Ci - ^ 



,(-5(V^'-16Ci+A)t) 



A2 - 16Ci 



(A15) 



This expression coincides with the one derived by 
Itakura and Tokura^ foUowing an elaborate Dyson series- 
type time-domain treatment. The off diagonal term 

initially follows a Gaussian decay, | {pti{z))^ {t)\ ^ 

i exp {— 2C^i2}, discussed in [7], and it becomes ex- 
ponential in the long time regime. The time constant 
appearing in the exponential decay is referred to as 
T2 following the NMR literature ^, turning out to be 
for Ca << A. Transforming to the dot 



basis the behaviour of the elements of the density opera- 
tor is exhibited in Fig. (A. 2). 




FIG. 2: {pniz))LLit) (Solid), {pniz))nRit) (Dashed), 
and [{pTi{z))^j^ (t)] (Dotted) are plotted as a function of 
time t for A = 1, = 1-5 and A = 0.5 
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